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Abstract 

Based on a canonical approach and functional-integration techniques, a series expansion of Green's 
function of a scalar field, in the presence of a medium, is obtained. A series expansion for Lifshitz- 
energy, in finite-temperature, in terms of the susceptibility of the medium is derived and the whole 
formalism is generalized to the case of electromagnetic field in the presence of some dielectrics. A 
covariant formulation of the problem is presented. 
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I. INTRODUCTION 



Quantum field theory is the quantum mechanics of continuous systems and fully developed 

j i 

in quantum electrodynamics which is the most successful theory in physics Quantum 
field theory trough path-integrals bridge to statistical mechanics and its applications include 
many branches of physics like, particle physics, condensed-matter physics, atomic physics, 
astrophysics and even economics Usually we are interested in a quantum field which 
has to be considered in the presence of a matter field described by some bosonic fields. For 
example, in quantum optics there are situations where the electromagnetic field quantiza- 
tion should be achieved in the presence of a general magnetodielectric medium 13-51 or in 
calculating the effect of matter fields on Casimir forces [6|, |7j. In these cases the matter 
field should be included directly into the process of quantization. Unfortunately there are 
very few problems where the interested physical quantities, like for example, the Casimir 
force, can be determined analytically and so finding an effective approximation method 
is necessary. A fundamental quantity in a quantum field theory is the propagator or the 
Green's function js| from which many physical quantities may be extracted. Here, using 
path-integrals, and based on a microscopic approach, we begin from a Lagrangian and ob- 
tain an expansion for the two-point correlation function i.e., the Green's function in terms 
of the susceptibility function of the medium for both scalar and electromagnetic fields in the 
presence of an arbitrary linear magnetodielectric medium. As an example of applications of 
these expansions for the case of a real scalar field we have introduced an expansion for the 
free energy or Lifshitz energy in the presence of some arbitrary dielectrics [9|, [l0| • Also we 
have considered the covariant formulation of the electromagnetic field in the presence of a 
linear magnetodielectric 



Casimir effects 



12|. 



llj which may have applications in quantum optics or dynamical 



II. SCALAR FIELD 



Let us start the section with a simple but efficient field theory which have a wide range of 
applications in many branches of physics, i.e., the Lagrangian of a real Klein-Gordon field 
in 3 + 1-dimensional space-time (x = (x, x°) £ M 3+1 ), with the following Lagrangian density 

C s = l -dMxW^{x) - \™\\x), (1) 
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and let the medium be modeled by a continuum of harmonic oscillators which is usually 
called the Hopfield model of a reservoir [l3] 

C m =\J™du (y 2 (x)-u 2 Y 2 (x)), (2) 

the interaction between the scalar field and its medium is assumed to be linear and described 
by 

poo 

£ in t= / duf{w,x)Y u (x)(p(x). (3) 
Jo 

Having the total Lagrangian we can quantize the system using path-integral techniques. An 
important quantity in any field theory is the generating functional from which n-point cor- 
relation functions can be obtained from successive functional derivatives. Here our purpose 
is to find two-point correlation functions or Green's functions in terms of the susceptibility 
of the medium. For this purpose let us first find the free generating functional which can be 
written as 



W [J,{J U }) = [ Dipeif d4x{ -^ [n+m2]v,+Jv} [ JJ J DF w es/^io 00 M-|^[^W]y.+ 



(4) 



where we have defined 



A=-(n + m 2 ), Bui = -(d 2 + oj 2 ), 

p(x) = -J(x), p u (x) = -J u (x). (5) 

Now using the following formula 

D<p(x)e-*M A M + W<p) = (detA)-h^ plA ' llp) (6) 

Eq. (j4]) can be rewritten as 

W [J, {J w }] = Ar e f<^>> e U~^WI^V> ( 7 ) 

where N = (detA)~^ W^idetBj) - ^ is a renormalization factor. Also from the following 
definitions 

Un + m 2 )G°(x,x') = 5\x-x') 
n 

Ud 2 +io 2 )Gl(x,x') = 5\x-x') (8) 



we will find 

j4l ik-(x-x') 

A' 1 = G°(x,x')=ih ' 



(2tt) 4 k 2 - m? 

2^ (k°) 2 -u 2 



B- 1 = G£(s, x>) = zM 3 (x - x') / (9) 



with the following Fourier transforms 



G°(k) lh 



k 2 — m? ' 

respectively. The free generating functional can now be written as 

W [J, {J w }] = jv e -iM d M dVj{x)G0{ ^ (11) 

and the interacting generating functional can be obtained from the free generating functional 
using the following formula [8[ 

W[J,{J U }] = e* /d4s/ "^ /( ^ ) (^)^(*^)w [J,{Ja,}] 



(12) 

Having the generating functional, the two-point function, i.e. the Green's function can be 
obtained as 
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(13) 

Now let us assume that the coupling function between the Klein-Gordon field and its medium 



is weak-one can also assume that the susceptibility of the medium is not far from vacuum- 
such that it can be considered as an expansion parameter which can be used to find a series 
solution for the Green's function or the correlation function. Using Eq. (ll2j) and after some 
straightforward calculations, we find the following expansion for the Green's function in 
frequency domain 

G(x-x» = G°(x-x»+ I d 3 z 1 G°(x-z 1 ,o;)[a; 2 x(a;,z 1 )]G (z 1 -x',w) + 

Jn 

d 3 z 1 d 3 z 2 G°(x - zi, u)[co 2 x(uj, z 1 )}G°(z 1 - z 2 , co)[u 2 x(^, z 2 )]G°(z 2 - x', cu) H 

(14) 




n Jn 



and it can be easily shown that it satisfies the Green's function equation which we will find 
in the next section. Note that since x) = 0, for x ^ Q, so we can rewrite the expansion 
(I14p in a more compact or matrix form as fallows 

G(u) = G°(uj) + G°(u) [co 2 x(u)] G°(u) + G°(u) [u 2 xM]G°M [oo 2 xH] G°(u) + ■■■ 
= G\u)[I-u 2 x{u)G°{u)]-\ (15) 

A. Equations of motion 

In this section we find the equations of motion for the fields and in particular we obtain a 
Langevin type equation for the scalar field 

d »7v7r^ = 7T^( n + m2 ^= dwf{u,x)Y„{x), (16) 
d{d^) dip J 

By solving Eq. (ll7j) and inserting it into Eq.(13), we find a Langevin equation for the Klein- 
Gordon field 

(□ + m 2 Mx) + ^ J dt' X (t-t',x)-^ V (x,t')=Z(x), (18) 

where x{ T -> x ) is the susceptibility function or the memory of the medium with the following 
Fourier transform 







The source field £(x) is defined by 

POO 

£(x)= / dwf(u,x)Y^(x) (20) 

JO 



where 

Yffat) = cos( W t)F w (x,0) + ^(x,0). (21) 

From Eqs. (120f2ip . we see that the source field depends on initial values of the reservoir fields 
so it can be considered as a noise field. The Green's function of Eq. (ll8j) satisfies 

(□ + m 2 )G(x-x / ,t-t') + -?- f ^"x(f-i" > x)^-G(x-x',t"-n = ^x-x'.t-t'). (22) 

dt dt" 
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In a homogeneous medium, where the memory function is position independent, Eq. (l22|) can 
be solved easily in reciprocal space 

1 1 



G(k,u) 



k 2 - u 2 + m 2 - u 2 x(u) k 2 - u 2 + m 2 - f du> £ f £'L 

1 (23) 



k 2 — u 2 e(u) + m 2 

where e{u) = 1 + x(u) can be considered as the dielectric function corresponding to the 
medium. From Eq. fT23l) it is clear that the Green's function in the presence of a homogeneous 
medium can be obtained from the Green's function of the free space simply by substituting 
uj 2 with uj 2 e{uj). Eq. ffTgj) in frequency-space can be written as 

(- V 2 - u 2 e(u, x) + m 2 ) G(x - x', u) = 5(x - x') (24) 

In some simple geometries the dielectric function e(w,x) depends on x as follows 

{e(ui) if x E £1 
(25) 
i if x ^ n 

where Q is a region or the union of regions where the space is filled with a homogeneous 
but frequency dependent medium with the dielectric function e(ou). In this case the Green's 
function can be found in some regular geometries [6] but for an arbitrary dielectric function 
it is quite complicated and in this case we find a series solution in terms of free Green's 
function and susceptibility of the medium. Note that in some geometries, electromagnetic 
field can be considered as two massless Klein-Gordon fields, and the scalar formalism can 
he.p for e Xam p, in obtaining LifsMtz energies or Ca 8 i m ir forces in snch geometries fl. 

B. Partition function 

Having the expansion (114"]) let us find the partition function in the presence of some dielectrics 
defined by the dielectric function e(u,x) which as a special case may be given by (125]) . The 
partition function of a real scalar field in the presence of a medium according to the modified 
Green's function given by ( |24l) can be written as 

~ = | D<peV = j D V e^ d4xC (26) 
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where C is given by (pQ). The partition function in frequency domain can be written as [ijj 

E = J D(fe~^f^ Jd 3 x^(x,-u;)[-a; 2 e(^x)-V 2 +m 2 ]^(x^)_ ^7) 

If we make a Wick rotation uj = %v in frequency domain the action will be Euclidean and 
the free energy can be determined from E = — ^ InS, where r is the duration of interaction 
which is taken to be sufficiently large. Using standard path-integral techniques we will find 
the free energy in finite temperature T as 

oo / 

E = k B T In det[^(^; x, x'] (28) 
1=0 

where V\ = Inlk^Tjfi is the Matsubara frequency, kg is the Boltzman constant and the 
prime over the summation means that the term corresponding to Z = 0, should be given 
a half weight. The kernel K(uf,r,r') = [ufe{ii>i, r) — V 2 ]5 3 (r — r'). Using the identity, 
In det[i^] = tr ln[K] and the fact that K(uf, r, r') = G~ x (ivi\ r, r') we find 

oo / 

E = -k B T Y ^ ln[G(iVi] r, r'] (29) 
z=o 

now using the expansion (fT^j) . we find the following expansion for free energy in terms of 
the susceptibility 

E = k B Tj2J2 / dhl • ' ' dhn ri ~ r2 ) ' ' ' G V*' r "- ri ) n) • • • Xfa, r«) 

1=0 n=l "* 

(30) 

where G°(r — r'; iv{) is given by 



i P ~\/ m2 + u ? I r - r 'l 

r - rO = T . (31) 

Air r — r 



which corresponds to a Yukawa potential with the modified mass ^m? + vf. 
III. ELECTROMAGNETIC FIELD 

In this section we use the Coulomb gauge i.e. V • A = 0, A = and find a similar expansion 
for the Green's function of the electromagnetic field in the presence of some arbitrary regions 
of matter which as an example can have applications in calculating the Casimir forces. For 



this purpose let us take the total Lagrangian density as follows 15] 

1 11 /*oo f 

C = -e (E 2 - — B 2 ) + - / dw{Y u (x) - co 2 Yj(x)) + / duf{u>, x) A ■ Y w (32) 
z Ho * Jo J 



The interacting generating functional can be written in terms of the vector potential and 
the medium fields as 



W = J D[A] Jj£)[y w ]expi J tfx^-^AikyA 



OO 1 

,, , -Y^(d 2 + u 2 )5 i3 Y u 
o 1 



+ / duf(oj,x)A i Y UJii + J i A i + / duJu,iY Wj ,, 
Jo Jo 

where summation over repeated indices is assumed and the kernel Ki~ is defined by 



(33) 



K, 



e d 2 - -V 2 



5ij H <9j<9j 



Now from the equation 



(34) 



(35) 



2 5Ji(x)5Jj(x') 

and similar calculations we will find the following expansion for the Green's function in 
frequency domain 

G fj (x -<£,;) = Gj(x-x» + / d 3 ziG5(x-zi,o;)[a; 2 x(tc;,zi)]Gj5(zi -x',w) + 

in 

/ / rf 3 zirf 3 z 2 Gi°(x- zi,w)[w 2 x(w,zi)]G z ° m (zi - z a ,w)[w 2 x(w,z a )]G^.(z a — x', to) H 



(36) 



which in matrix form can be written as 



G(w) = G°(w) + G°H[w 2 x(^)]G°(a;) + G°M[w 2 xM]G°(w)[w 2 xM]G°M + • • • 

= G°M [l - uj 2 x(u) G°(w)] _1 . (37) 

A. Equations of motion 



From Lagrangian density (j2J) we find the following equations 



KijAj 



Y u>i + w 2 Y w ,i = -f(u 1 x)A i 
Solving Eq. (!39|) and inserting the solution into Eq. (l38l we find 

KijAj + j t Tdu f 2 {u, x) / dt' G u {t - t')-^A(t>) = I du f(u, x)Y* 



df 



(38) 
(39) 

(40) 



where 



duj' e ^'(*-*') 



(41) 



2n uo 2 — cu' 2 

and Y^i is the homogeneous solution of Eq.( l39l) which depends on the initial values of the 
medium fields and can be considered as a noise or fluctuating field which does not affect the 
Green's function. Using Eqs.f l20p and ( T4T1) we can rewrite Eq.(28) as 

f) f°° r\( / C f) f°° 

KijAj + | / o ^ x) J dt> e^^^Mt') = jf du /( W> x)Y* (42) 

which in frequency- domain can be written as 



-e u 



1 -V 2 )^- • -<),<), 



/'() 



1 

/'o 



Aj(x ) u)-wx(w ) x)i i (w,x)= / dwza;/(w,x)l^.(w,x) 

'o 

(43) 



The Green's function of Eq. (143]) satisfies 



-e (l + x) w 2 Sij) - — V 2 ) <y y + —didj 



/'o 



/'o 



G^(x, x', u) = — <5 3 (x - x')<% (44) 



1*0 



which using the definitions e(w,x) = e [l + x(w,x] and c 2 = e /io can be written as 



to 



-o-e(o;, x)5ij - V % + <9i<9j 



G 3 -fc(x, x , w) = 5 3 (x - x')S, 



ik • 



(45) 



and it can be easily shown that the Green's function fl36l) satisfies Eq.( l4"5]) . 

A similar approach can be followed to find the partition function in terms of the suscep- 
tibility of the medium as follows 



1=0 n=l 



(46) 



where the free Green's function G?-(x — x'; iv{) satisfies Eq.f l45|) with e(u, x) = 1 and u = ivi. 
By defining r = x — x', we find 



J cr Attt 



5^(1 + — + 



1 2 ' 



TiTj .„ 3c 3c 

1 + — + 



2 2 • 



+ -^-5 3 (r). (47) 



In zero temperature, the summation over the positive integer I is replaced by an integral 
according to the rule h / — f{ iy ) ~~ ^ &bT f{^i)- For a nice discussion of Casimir- 

JU 1=0 

Lifshitz interaction between dielectrics of arbitrary geometry see 10]. 



9 



B. Covariant formulation 



In reference [ll|] electromagnetic field quantization in a moving medium has been investigated 
by considering the medium to be modeled by a continuum of tensor fields. But the medium 
can also be modeled by a continuum of scalar fields, i.e., Klein-Gordon fields which we will 
follow here. So we consider the electromagnetic field interacting with a moving medium, 
a situation which can have applications in dynamic Casimir effects. For this purpose we 
consider the following Lorentz invariant Lagrangian density 

-1 1 /*oo /»oo 

C(x) = -d li A v d*A" + -j dcu[d^Y u d^-co 2 Y 2 ] + J <hjr{u,*)Y^A v . (48) 

where / M "(a;,x) is an antisymmetric coupling tensor which couples the electromagnetic field 
to its medium and is related to the susceptibility of the medium trough Eq. ([52]) . From 
Euler-Lagrange equations we find 

DA v = -d li K' a '(x), (49) 
(□ + u 2 ) Y u (x) = r(oo, x)d ll A v {x) (50) 

/■oo 

where □ = d 2 — V 2 , and the antisymmetric tensor K^ v (x) = / du / M "(u;, x)Y u can be 

Jo 

considered as the polarization tensor of the medium. By solving Eq. (150j) and inserting it 
into Eq. (09]) we find 

DA^x) - J d^x' d u d' a x v ^{x,x')A p (x') = -d n K N ™ (51) 

where K N '^ U = J °° dui f^ u (co, x)Y^ , and Y^ is the homogeneous solution of Eq.( l50l which 
can be considered as a noise field. The susceptibility tensor x" ( x i x ') is defined by 

/"OO 

X v " af> (x,x')= duf^(u,x)Gl(x-x')f aP (u,x') (52) 
Jo 

where G° (x — x') is given by Eq.flU]). The Green's function in this case satisfies the following 
equation 

□ G^(x - x') - J d 4 x" g^ 5 d, d'i x j5aP (x, x") G^x" - x) = g, u 6\x - x') (53) 

and for the Green's function we find the following expansion in terms of the susceptibility 
tensor 

G^x, x') = G%(x d% d 4 z 2 GlJx - Zl ) r™(z u z 2 ) G° u Jz 2 - x') + 

J d 4 Zl ■ ■ ■ d 4 z A Gl Ul (x - z x ) T^{ Zl , z 2 ) G° U2U3 (z 2 - z 3 ) T^(z 3 , z 4 ) G° U4U (z, ~ x') + ■ ■ ■ 

10 



where for simplicity we have defined T UlU2 (zi, z 2 ) = d^d^x^^^^izi, z 2 ). So given the 
susceptibility tensor of the medium, one can find the Green's function perturbatively in 
terms of the susceptibility. Having the Green's function-at least perturbatively- we can 
investigate for example the dynamical energy configurations, which is closely related to the 
problems of dynamical Casimir effect, which is under consideration. 

IV. CONCLUSION 

Based on a Lagrangian approach, scalar and vector field theory in the presence of a medium, 
modeled by a continuum of Klein-Gordon fields, considered and a series expansion for the 
Green's function of the theory in terms of the susceptibility of the medium obtained. From 
the partition function of the scalar field, an expression for the free energy in terms of the 
susceptibility of the medium obtained and the formalism generalized to the case of elec- 
tromagnetic field in the presence of some dielectrics. Also, the covariant form of the elec- 
tromagnetic field in the presence of moving media investigated and an expression for the 
Green's function in terms of the susceptibility tensor obtained which can have applications 
in dynamical Casimir effects. 
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